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m 

During the past decade there have been several attempts to detect cosmogenic ultra 
high energy (UHE) neutrinos by searching for radio Cerenkov bursts resulting from 
charged impact showers in terrestrial ice or the lunar regolith. So far these radio searches 

i-C have yielded no detections, but the inferred flux upper limits have started to constrain 

physical models for UHE neutrino generation. For searches which use the Moon as a 
target, we summarize the physics of the interaction, properties of the resulting Cerenkov 
radio pulse, detection statistics, effective aperture scaling laws, and derivation of upper 
limits for isotropic and point source models. We report on initial results from the 
RESUN search, which uses the Expanded Very Large Array configured in multiple 
>• sub-arrays of four antennas at 1.45 GHz pointing along the lunar limb. We detected 

no pulses of lunar origin during 45 observing hours. This implies upper limits to the 
differential neutrino flux E 2 dN/dE < 0.003 EeV kuT 2 s" 1 sr" 1 and <0.0003 EeV km~ 2 
s _1 at 90% confidence level for isotropic and sampled point sources respectively, in the 
neutrino energy range 10 21 ' 6 < E{eV) < 10 22 ' 6 . The isotropic flux limit is comparable 

On to the lowest published upper limits for lunar searches. The full RESUN search, with an 

additional 200 hours observing time and an improved data acquisition scheme, will be 
be an order of magnitude more sensitive in the energy range 10 21 < E(eV) < 10 22 than 

^ previous lunar-target searches, and will test Z burst models of neutrino generation. 



Subject headings: ultra-high energy neutrinos, cosmic rays, lunar interactions 



1. Introduction 

Determining the origin and characteristics of extragalactic cosmic rays is of fundamental im- 
portance in astrophysics. Direct detection of high energy cosmic rays from distant sources is difficult 
due to interaction with cosmic microwave background photons. Above a limiting energy of 10 19 ' 5 
eV, cosmic rays interact with CMB photons within a few tens of Mpc, producing a high-energy 
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pion shower (Greisen 19661. The pions quickly decay, producing ultra-high energy (UHE, E > 
1 EeV = 10 18 eV) neutrinos that can propagate unimpeded directly to the Earth. This limit, 
known as the Greisen- Zatsepin-Kuzmin (GZK) limit, predicts that UHE particle processes thought 
to occur in distant AGN cores and hypernovae cannot be seen directly at Earth, but must be traced 
via the secondary production of UHE neutrinos. Indeed, recent measurements of the cosmic ray 



spectrum at UHE energies (Yamamoto 20081 show a clear drop above the GZK limit, consistent 



with the GZK prediction. Other, more speculative, sources for UHE neutrinos include decaying 
super-massive particles, such as magnetic monopoles, and cosmic topological defects ( Stanev|2004 l. 
Although not yet detected, these UHE neutrinos will provide an fundamentally new window for 
studying the physics of AGN and possibly exotic astrophysical particles. 
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Fig. 1. — Observed Cosmic Ray spectrum across a wide range of energies dCronin et al.|1997[ ). Inset 
shows a detailed view of the "Ankle" region located at the Ultra High Energy end of the spectrum, 



from Waxman ( 2009 1 . Recent measurements at UHE energies have revealed a drastic attenuation 
of the CR spectrum at IO 19 6 eV, consistent with the predicted GZK cutoff. The ordinate of the 
inset has been scaled by a factor of E 2 ^ 1 to better illustrate the CR flux drop-off. 
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Since the interaction cross-section of neutrinos with matter is notoriously small, large targets 
are required to detect UHE neutrinos. A promising approach is the observation of radio bursts 
resulting from UHE neutrino interactions with the Moon (Fig. [2]). Askaryan (19621 calculated that 
UHE particles would interact with bulk materials such as ice, salts, or the lunar regolith. The 
regolith interaction produces a shower of particles with a resulting a 20% electron excess. As the 
electrons are moving at the free-space speed of the incident neutrino ~ c, larger than the phase speed 
of light in the medium, they produce a very short duration (hence wide bandwidth) burst of radio 
Cerenkov radiation. The radiation is concentrated in a thin cone visible by cm- wavelength radio 
telescopes. Askar'yan's prediction and many details of the particle interaction have been directly 
verified with pulsed electron bunches and salt and ice targets at the Stanford Linear Accelerator 
( Saltzberg et~aT|2001| piocinovic et aT|2006| [Gorham et al|2007| ) . The pulse emission was shown 
to be broadband, linearly polarized, with nanosecond characteristic timescales, and corresponding 
bandwidth 1-10 GHz, depending on the observer's location in the radiation cone. 



Interaction w/ 




Fig. 2. — UHE Neutrino interaction with lunar regolith matter. A neutrino-parton interaction 
results in a hadronic shower. The shower creates an energy cascade of e + e~ pairs with a 20% charge 
excess and consequent Cerenkov radiation cone. The cone partially escapes the lunar surface, where 
an observer sees a nanosecond-duration radio pulse. 



Searching for radio bursts with such small time scales is problematic. Although the peak flux 
density is predicted to be very intense, e.g. of order 10 5 Jy at 1 GHz ( Alvarez-Muniz et al. 20001, 
interference and radiometer noise statistical fluctuations can mimic the desired signals. Hence, 
well-designed detection schemes must strongly discriminate signals of lunar origin from 'accidental' 



(statistical fluctuation) pulses. The first pioneering search was performed in 1995 when Hankins 



et al.| (1996) observed the Moon in a band from 1200 -1700 MHz using the Parkes telescope in 
Australia. They searched for pulses dispersed by the Earths ionosphere, using the ionospheric 
dispersion delay to discriminate against terrestrial RFI pulses. A large fraction of their 12 hour 
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experiment was spent pointed at the lunar center, not realizing that only on the lunar limb could 
Cerenkov emission escape in the direction of the observer. 

The second lunar radio search, known as GLUE (Goldstone Lunar UHE Neutrino Experiment, 



Gorham et al. (2004)) employed the Goldstone 70m and 34m telescopes with a large (22 km) 
separation to search for delayed pulse coincidence. They also found no pulses of lunar origin in 
120 total hours of observation, but the sampled fraction of the lunar limb was limited by the small 



primary beams of the antennas compared with the angular size of the Moon. Beresnyak et al. (2005) 



described a search using the 64-m diameter Kalyazin telescope using four simultaneous frequency 
bands between 0.6 and 2.3 GHz with the telescope pointed along the lunar limb. They observed 
for a total of 31 hours, but detected no pulses with the expected frequency-dependent differential 



delay from a lunar source. Scholten et al. (2009) reported on first results from the 'nuMoon' search 



using the Westerbork Synthesis Radio Telescope in several bands within the frequency range 115- 
180 MHz. They had no detections in a preliminary analysis of the first 20 hours of observation, but 
are continuing to observe. Future lunar radio searches using existing or planned radio telescopes 



have been discussed by Scholten et al. (2006 LOFAR), Panda et al. (2007 GMRT), |James k 



Protheroe (2009 SKA), and James et al. (2009, ATCA). Finally, we note that several terrestrial 



searches for Cerenkov radio bursts expected from neutrino interactions with the Earth's polar ice 



sheets have been reported e.g. FORTE (Lehtinen et al. 2004), RICE (Kravchenko 2006), ANITA- 



lite ( |Barwick et "aL]|2006| ) and ANITA ( |Gorham et~aI"1|2009D . So far, no cosmogenic UHE neutrinos 
have been detected by any of these experiments. 

In this paper we describe the physical and technical considerations relevant to searches for 
Cerenkov radio pulses resulting from UHE neutrino interactions with the Moon. Although several 
previous papers have described various aspects of this scheme, our goal is to fully describe the 
relevant physical and technical calculations in a self-contained manner. In addition, we apply 



scaling laws derived from the recent analytic effective aperture calculation of Gayley et al. (2009) 



and introduce a new graphical representation of non-detection upper limits which better accounts 
for the non-detection probability as a function of predicted model flux. 

We first discuss the neutrino- regolith interaction physics, properties of the resulting Cerenkov 
radiation, the effective detection aperture, and statistical considerations for multiple antenna de- 
tection schemes. We apply these to derive the sampled neutrino energy range and isotropic flux 
as a function of the energy-dependent pulse profile and detector parameters. We then describe 
the RESUN search (Radio EVLA search for UHE Neutrinos), a radio search using the EVLA0 
We summarize the experimental setup, including the use of multiple sub-arrays and high-speed 
pulse detection hardware to optimize sensitivity to short-timescale pulses. Finally, we report initial 
results of the search, and compare them with previous radio searches. 



1 Expanded Very Large Array, near Socorro NM, operated by the National Radio Astronomy Observatory. 
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2. Physics of the UHE Neutrino Lunar Interactions 
2.1. Neutrino-Regolith Interaction 



Incident UHE Neutrinos interact with baryons in the lunar regolith with an energy-dependent 

E 



cross-section (Reno 2005) 



a{E) = 1.57 x 10 



-31 



ZeV 



1/3 



cm 2 . 



(2.1) 



where ZeV = 10 eV. The corresponding mean free path of the neutrino can be expressed as 



a{E)p 



/ZeV\ 1/3 
35.7km f ——J 



(2.2) 



where E is the neutrino energy and p = 1.8 gm cm 3 is the mass density of baryons in the regolith 



(James & Protheroe 20091. The neutrino does not impact baryons per se, but instead with their 



interior, point-like constituent partons (quarks and/or gluons). The interaction results in a hadronic 
shower with approximately 20% of the initial neutrino energy. This shower quickly evolves into a 



short-lived current source via electron entrainment and positron annihilation (Alvarez-Muniz et al. 



2008). This transient current creates a pulse of coherent Cerenkov radiation, since the electrons' 



speed (c) is greater than the local phase speed of light in the medium. 
The pulse is emitted in a cone of half-angle 

6 C = acos( — ) ~ 55° 
n r 



where n r = 1.73 is the refractive index of lunar regolith (James & Protheroe 2009). The pulse is 
broadband with a characteristic timescale r = c/L s , where the shower length L s in regolith is given 



by (Scholten et al. 2006) 



UE) - 1.-0 1 +0.08 log ( m 



(2.3) 



The corresponding shower timescale is r ~ 5 nsec with a small dependence on energy. The remain- 
ing 80% of the neutrino energy goes into one of two different reactions. For neutral-current (NC, 
branching ratio ~ 1/3) interactions, the neutrino is scattered and retains the remaining energy. In 
the charge-current (CC, branching ratio ~ 2/3) case, a muon, electron, or tau lepton is produced 
with the remaining ~ 80% of the energy. While muons and tauons are absorbed by the surround- 
ing medium, electrons will initiate a electromagnetic cascade which in turn creates a secondary 
Cerenkov cone. However, the EM cascade is considerably elongated by the Landau-Pomeranchuk- 



Migdal (LPM) effect (Alvarez-Muniz et al. 2000), forming a narrow emission cone which has a 
negligible contribution to the total emission at most observing angles. 
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Fig. 3. — Sketch of the neutrino-parton interaction. The hadronic shower carries 20% of the 
neutrino energy, eventually in e + e~ pairs. The positrons are lost more quickly, resulting in a net 
20% e~ charge excess. The charged current emits a Cerenkov pulse with a cone width given by 
equation |2.4 



2.2. Properties of the Cerenkov Radiation Pulse 

The hadronic shower creates a Cerenkov radiation cone whose peak amplitude and angular 
width depend on energy E and observation frequency v. The cone's full angular width at e _1 in 
lunar regolith can be written (James &; Protheroe 2009) 



A0 H (E,u) = 2.9° 



"GHz" 




V 





1.15 + 0.075 log 



(0.2E 

\ZeV 



(2.4) 



where we have multiplied the James & Protheroe (2009) expression by 1.2 to convert from half- 
wifth to e _1 width. At E = 1 ZeV, the resulting cone widths vary from 17.6° at 150 MHz to 1.3° 
at 2 GHz. 

In the absence of regolith attenuation, the peak electric field strength {9 = 9 C ) of a lunar 
Cerenkov pulse received at the Earth can be parameterized as (Gorham et al. 2000| James &z 



Protheroe 2009) 



^(^,^ = 0.51 



0.2E 

ZeV 



v 



1 + 



m" 1 MHz 



(2.5) 



where vq = 2.32 GHz, a = 1.23 for regolith material (James &; Protheroe 2009) and v is the 
observation frequency. For Phase-A of the RESUN search (fo=l-45 GHz, Av = 50 MHz), the peak 
electric field of the Cerenkov pulse at the telescope, including transmission loss the lunar surface, 
is 



S c max {E) = 0.041 



\Zev) 



HV m- 1 MHz" 1 



(2.6) 
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and the e 1 width (ignoring the small effect of the sin terms in 2.9) of the Cerenkov cone is 
A6 H ~ 2.1° for all E > E . 

Unfortunately, the burst radiation will be strongly attenuated by the local medium as it prop- 
agates to the surface from the shower interaction in the interior. The e _1 attenuation length for 



the electric field in lunar regolith is (Olhoeft & Strangway 1975) 



X e = 18 m • 



GHz 



(2.7) 



Hence, the electric field strength from each escaping ray path along the Cerenkov cone must be 
corrected by a factor 

rz sec(9) 

£{z) = £(0) / e~ {s/Xe) ds , (2.8) 
Jo 

where z is the vertical depth of the shower and 9 is the ray path angle to the surface. The 
short attenuation length implies only glancing angle impacts produce observable Cerenkov cones. 
Neutrinos with large downward angles create radiation cones which are completely absorbed or 
suffer total internal reflection. 



The angular dependence of the Cerenkov cone electric field is ( Alvarez-Muhiz et al. 1 |2006| 



where 



S c (6,E,u) = S c max (E,u) e -C(^) 3 



C(9,E,u) 



sin(0) 



sin 



(2.9) 



A8 H (E,v) 



Equation 2.9 represents the burst electric field profile interior to the lunar surface. Of greater 
interest is the transmitted electric field after propagating through the regolith-vacuum interface. 
The transmission across this interface is a complicated function of both incident angle and polariza- 
tion (Williams 2004 Gusev et al.||2006 ). However, Gay ley et al. ( 2009 1 argue that the transmission 
can be approximated by a constant value tu ~ 0.6 without significant error. With this in mind, we 
write the maximum free-space electric field strength as 



£f(E,v) = 0.6£ c max (E,v) 



(2.10) 



2.3. Signal-to-Noise Ratio and Minimum Detectable E-Field 

must exceed the minimum 



2.5 



To detect the Cerenkov burst, the electric field given by Equation 
receiver sensitivity for each detector element (single antenna). Note that this is different from 
interferometer arrays, in which signals are co-added coherently. The RMS electric field fluctuations 
at the feed of a radio telescope with effective area A e , system temperature T sys , and bandwidth Au 
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can be calculated by equating the power associated with a radiation field having and RMS electric 
field £ illuminating the telescope area A e with the shot noise power of the receiver, 

— Au = — 

Z A e 

so that the RMS field due to detector noise fluctuations is given by 

^(„W) = ( , n) 

where rj is a dimensionless constant which depends on the polarization properties of the telescope 
feed and radiation, kb is Boltzmann's constant, and Zq = 377 f2 is the impedance of free space. For 
a given detector scheme, the minimum detectable electric field is a multiple N a of the RMS field 
given by the detection threshold as described in section |2.4| For Phase-A of the RESUN search, 
(N& = 3.98,7] = 2 for linearly polarized radiation and circularly polarized feeds, T S y S ~ 120 K 
pointing on the lunar limb, A e = 343 m 2 ) this evaluates to 

= N a -£ d = 0.034 fiV nT 1 MHz" 1 (2.12) 



By equating equations 2.10 and 2.12 we find that the minimum detectable incident neutrino energy 



for RESUN-A is E mir> = 1.38 ZeV. 



2.4. Setting the Threshold: Accidental Trigger Rates 

At each detector, a threshold level must be set to discriminate against false pulse detection 
arising from statistical noise fluctuations. Assuming the receiver noise voltage obeys Gaussian 
statistics and that the receiver maintains a linear response within the voltage range of the fluctua- 
tions (see section 4) , the probability of a sample signal a exceeding a specified threshold at is given 
by: 



p Q {a > a t ) = -j= J exp\-—jda = ei{c\-j=) (2.13) 



where erfc is the complementary error function. With fast sample times, the false signal rate can 
be quite high, even for very high threshold values. For example, with a 100 MHz sample rate (10 
ns samples) and a 6a threshold, one expects an 'accidental' trigger approximately once every 5 
seconds. 

The primary beam of a radio telescope will illuminate an area on the lunar surface with a 
range of geometrical delays corresponding to a window of N s time samples. Assuming po « 1, 
The probability p w of a sampled time exceeding the threshold level in the window of N s samples is 
simply 

p w {cj t , N„) = l-[1- Po(<Tt)] Ns ~ N s ■ Po {a t ) (2.14) 
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For an array of N a telescopes whose primary beams are coincident on the Moon, the joint probability 
Pj that a statistical fluctuation will exceed the threshold level at all telescopes is the product of 
the individual antenna detection probabilities, where in general, each antenna may have a different 
signal to noise ratio and corresponding threshold level. Thus, the joint probability of an accidental 
trigger in time window of N s samples (after correction for geometrical delays) for N a antennas is: 

N a 

Pj (a t ,N s ,N a ) = A^" 1 ) l\po(au) (2.15) 

i=l 

where we have accounted for the possibility of different threshold levels an for each antenna. 
The result is a dramatic decrease in the accidental trigger rate, even with small arrays of a few 
telescopes. For a given experiment, the threshold level must be set high enough that the probability 
of an accidental trigger during the observation is highly unlikely. 




Fig. 4. — Expected number of accidental detections per second versus threshold level calculated 



using equation 2.15 Accidental rates are shown for arrays of one to six telescopes using 10 ns 
sampling time and 140 ns. delay tolerance window. The measured accidental detection rate for five 
sample one-minute intervals (scaled to 50 hours) for one telescope (x's), two telescopes (squares), 
three telescopes (triangles), and four telescope (upper limit) detection schemes are shown using a 
4o" threshold are in excellent agreement with the noise model. 
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2.5. Expected detection rate: Effective Aperture Calculation 



Calculation of the expected detection rate is a difficult problem because of the non-trivial 
geometry of the partially escaping radiation cone, attenuation in the lunar regolith, the effects 
of random surface roughness, the frequency-dependent detector sensitivity, and the 2-dimensional 
dependence on incident neutrino angle. Several groups have solved this problem for specific ex- 
periments by developing Monte Carlo simulations, following ray paths for an variety of neutrino 



flux models ( 


Gorham et al. 


2000 


Protheroe 


2009 


). More recently, 



More recently, we have developed an analytic approach that expresses the ex- 
pected detection rate as a function of characteristic physical parameters ( Gay ley et al.|2009 1 , which 
we summarize here. 

The complete system of lunar target and detector can be parameterized by a direction- 
dependent detection aperture, which when multiplied by the incident neutrino flux, gives the de- 
tection rate in that energy bin and direction. In the following, we consider only isotropic neutrino 



fluxes, although this approach can be generalized to a direction-dependent case (e.g., James et al 



2007). 



The effective aperture A e can be written as the product of the cross-sectional area of the 
Moon, the total solid angle of incident neutrinos (4-7T for isotropic neutrino flux), and a detection 
probability P which takes into consideration the interaction specifics, 



A e (E,v,£i in )=A Q -P(E,v,£° 



(2.16) 



where E is the neutrino energy, v is the observing frequency, is the minimum detectable 



electric field strength of the detector (see Section 2.3), and Aq is the geometric lunar aperture 
QWmiams||2004l ), 



A = 4tt ■ (ttR 



2 ) 

ml i 



where Rm is the radius of the Moon. Note that the aperture has units of physical area times solid 
angle. The probability P is the fraction of neutrinos entering the Moon at energy E whose radio 
pulse will be detectable with a telescope (detector) having minimum electric field sensitivity £, 



d 

rain 



at frequency v. We now summarize the calculation of P detailed in Gayley et al. (2009). 
An approximate expression for P can be written 

16 16 



P = P ( f O A0 H + 



3vr 3 / 2 



cr s + 



-«0 



(2.17) 



where 



nr. 



1 



8n r 



and the three bracketed terms represent the contributions from downward incident neutrinos on a 
smooth surface, downward neutrinos on a rough surface, and upward traveling neutrinos respec- 
tively. In the first term, the dimensionless scaling parameter /q is the ratio of the thickness of the 
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Cerenkov cone at the e 1 full thickness to the thickness at the minimum detectable electric field 



e d min (cf. Fig. 2, Gayley et al. (20091, 



I hi 



(2.18) 



In the second term, a s is the Gaussian half-width RMS surface roughness angle of the lunar surface. 
A self-similar (fractal) model provides a good fit to radar surface return measurements, which we 



can write as a function of observing frequency (Shepard et al.||1995 Gayley et al. 20091 



cr (A) = V2 tan _1 (0.14 u - 22 ) ~ 0.2 v 



0.22 



(2.19) 



where a Q is in radians and v is in GHz. For the RESUN search {y = 1.45 GHz), we obtain 
a s = 12.2°. In the third term, ao is the maximum acceptance angle for upward incident neutrinos 
which contribute to the escaping radiation cone, and is given by 



a (E) 



IAX U (E) 
2R M 



1.8° 



E 



(2.20) 




Fig. 5. — Aperture versus detector sensitivity and observing frequency (v) for a fixed 

neutrino energy E = 10 21 ' 5 eV. Values for experiments RESUN Phase-A , RESUN Phase-B , 
GLUE, nuMoon, and SKA are also shown. 



-12- 




1Q 19 1Q 20 10 21 1Q 22 1Q 23 1Q 24 

Neutrino Energy (eV) 



Fig. 6. — Effective aperture versus neutrino energy for fixed detector sensitivity E m i n = 0.01 
/j,V m _1 MHz~ 1 and observing frequencies between 0.1 GHz and 10 GHz. 



Fig. [HJshows the effective aperture as a function of neutrino energy for fixed detector sensitivity 
(E min = 0.01 \N m _1 MHZ" 1 ) and a range of observing frequencies from 0.1 GHz to 10 GHz 
using equation 2.17| Note that the sharp low-frequency cut-off is a result of /o — > , which for 
fixed telescope sensitivity occurs at higher neutrino energy as the observing frequency is lowered. 
Frequencies near 1 GHz provide the best compromise between low-energy cutoff and maximum 
aperture above the cutoff. 



2.6. Flux Limits for Non-Detection 



The expected number of detections N v in exposure time t is the product of the neutrino flux 
and the effective target area integrated over the total observing time, solid angle, and all detectable 
energies, 

ft r poo 

N v (t,n)= / dt' dfl dE I(E,h) A(E,n(t')) , (2.21) 

JO J 4tt J E min 

where A is (units of area) is related to the aperture A e (units of area-sr) by 

A e (E) = I A(E, h) dn (2.22) 
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and I(E, h) is the differential energy neutrino flux in direction h. For the case of an isotropic 
neutrino flux, we have 

I(E,h)=I(E) 

where 1(E) is the flux per unit steradian. This results in an expected count in time t 

rt roc p 

N lso = dt' dE 1(E) / dn A(E, h) , 



I dt' / 

" E m i n 



■lir 



or, using eqn. 



2.22 



poo 

N iso = t dE 1(E) A e (E) . 
Alternatively, for a point source the neutrino flux is a delta- function in direction, 

I(E, h) = 1(E) 6(h-h s ) 



(2.23) 



(2.24) 



where h s is the unit normal direction to the source. Substituting into equation 2.21 we have 



I p = ['dt' f 

JO J E, 



V„ - / dt' I dE 1(E) A(E,h s (t') ^t p dE 1(E) A(E) , (2.25) 

J E m i n J E m i n 

where t p is the total time that the aperture intercepts the source as the Moon moves across the 
source region. 

In order to proceed, we need to consider the energy dependence of the aperture solid angle. The 
effective aperture directional distribution on the sky is a complicated calculation which has only 



been done using Monte Carlo simulations (James et al. 2007). However, we can use the analytic 



results of Gayley et al. (2009) to approximate the solid angle (SIa) by noting that at high frequen- 
cies iy > 1 GHz) the downward neutrino acceptance angular spread is dominated by the surface 



roughness angle <To (eqn. |2.17 ) and is much larger than the upward-directed contribution. There- 
fore the radial angular spread is approximately the surface roughness angle, which is independent 
of neutrino energy. We then approximate A as 



A(E) 



ME) 



where is an equivalent-width estimate of the aperture solid angle, which can be approximated 
as 



Collecting terms, we have the expected count rate for a point source 

f oo 



AT ~ llL 



poo 

/ dE 1(E) ME) , 



(2.26) 



Once the energy-dependent aperture A e (E) is known, equations (2.23) or (2.26) can be dif- 
ferentiated and used to determine the differential neutrino flux (dN/dE) as a function of energy. 



- 14 - 



For experiments with a null detection, at 90% confidence level the neutrino flux upper limit is 2.3 
times the differential flux upper limit, assuming Poisson counting statistics. Using to the commonly 
plotted quantity F(E) = E 2 dI(E)/dE, the upper limit for a non-detection in a given observing 
time t is 



F iso (E,t) < 2.3 



1 E 



t A e (E) 

while a single point source of neutrinos, the corresponding upper limit is 

n A E 



F p (E) < 2.3 



t P A e (E) 



(2.27) 



(2.28) 



where t p is the total time that the source is within the aperture solid angle. For lunar search exper- 
iments, the beam tracks the Moon's motion on the celestial sphere, so that the mean observation 
time can be written 



tr. 



105" 



where 6t (radians) is the effective angular extent of the aperture solid angle as it transects the point 
source. If the aperture is dominated by the lunar roughness term in equation 2.17, we can write 



do . 



Collecting terms and expressing the frequency dependence explicitly, we can recast eqn. (2.28) 
using approximate expressions above for t p , a Ql and 6 eq to obtain the flux detection limit (90% 
confidence) for point sources 



F p {E,u) ~ 1.2 x 10" 



GHz 



0.22 E 



A e (E) 



(2.29) 



3. RESUN Search 

The RESUN search uses the National Radio Astronomy Observatory's Very Large Array (VLA) 
to search for Cerenkov radio bursts originating from the Moon. The search comprises two phases: 
RESUN-A, a 45-hour search completed in February 2008, and RESUN-B, a 200-hour search with 
significantly enhanced limb coverage, wider total bandwidth, and lower pulse detection threshold. 
This section describes the experimental design and effective apertures for both searches and reports 
on results for the RESUN-A search. RESUN-B is scheduled for completion in early 2010. 



3.1. Array and receiver configuration 

The VLA consists of 27 25-m diameter parabolic reflectors with Cassegrain optics in an array 
whose maximum baseline varies from 1 km to 30 km depending on array configuration. For RESUN- 
A we used one sub-array of four antennas with antenna separations ranging from 0.6 km to 5.3 
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km. We observed right-circularly polarized radiation at 1.45 GHz center frequency with a 50 MHz 
bandwidth. Each antenna in the subarray was pointed to a target on the east limb of the Moon 
(FigJT^,). The VLA antenna FWHM primary beam at 1.45 GHz is 30 arcmin. Approximately 
one third of the lunar limb was illuminated in this configuration. Note that each antenna is an 
independent pulse detector so the array is not being used in the usual interferometer mode. 



Phase A: 

33% Limb Coverage 




Phase B: 

ioo% Limb Coverage 



X4 




|- 1 | 






Beams 



(a) 




Fig. 7. — Subarray beam geometries for RESUN-A (left) and RESUN-B (right) experiments. 



The VLA receivers and correlator are currently being upgraded to improve sensitivity and 



frequency agility (Butler et al. 2006). The RESUN-A search was scheduled during the transition 
from the 'classic' VLA system to the upgraded EVLA system. In this hybrid transition mode, both 
VLA and EVLA antennas were available. In order to determine the most appropriate system for 
RESUN, we measured the receiver output voltage dynamic range for both VLA and EVLA receiver 
systems by comparing histograms of sampled receiver output voltage to a Gaussian statical model 
(Fig. wj). A large voltage dynamic range is critical to Cerenkov pulse detection experiments, since 
the calculated pulse amplitudes much larger than typical astronomical sources. The exact threshold 



value depends on the experimental setup, as discussed in section |2.4| We found that the 'classic' 
VLA receivers were unusable for the RESUN search since they saturated well below the RESUN 
threshold of 4a, while the new EVLA receiver systems were in excellent agreement with Gaussian 
statistics until voltage excursions approached 5a. Hence we used a subarray comprised entirely of 
EVLA antennas. 



3.2. Pulse detection scheme 



The down-converted waveform at each antenna was sampled and processed using FPGA-based 



digital signal processors developed by the CASPER laboratory at U. C. Berkeley (Parsons et al 



2008). Each subsystem module consists of two dual channel 8- bit analog to digital converters 



(ADCs) attached to a Internet Break-Out Board (iBoB), which provides Internet protocol (IP) 
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Fig. 8. — Comparison of measured count rates for sample VLA (antenna 4, dashed line), EVLA 
(antenna 19, dashed-dot line) receivers, and Gaussian statistical model (solid line) versus threshold 
level for sampled voltage output of square-law detector (proportional to input power). Note that 
while the VLA receiver shows significant saturation above 3a, there is very little saturation for the 
ELVA system for signals < 5a, well above the RESUN threshold. 



based communication (UDP packets) with the host computer and allows user-configurable access 
to the connected devices. 

The raw input data stream consisted of four channels, one from each telescope, that were 
simultaneously sampled every 10 nanoseconds using a common sampling clock trigger. In order to 
avoiding recording all bits, we programmed the FPGA to only transfer samples (and time stamps) 
when the signal level in any one channel exceeded a predefined threshold level. A wave front 
originating from the Moon would illuminate all four antennas with differential time delays dictated 
by instrumental and time-dependent geometrical delays. Hence, we would record four threshold 
events per pulse, suitably displaced in time. A block diagram of the data acquisition system in 
shown in Fig. [9] 

We searched for 4-antenna pulse coincidences by calculating the expected geometrical delays at 
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Fig. 9. — Block diagram of RESUN pulse detection scheme for a 4-antenna sub-array. 



triggered event times in a single channel. We then searched for triggered events in the other three 
channels using delay windows centered at the expected total (instrument plus geometrical) delays. 
We calculated the geometric delay for sources along the illuminated chord of the lunar limb using 



the Chapront ELP2000/82 algorithm (Meeus 1998) which has < 10 arcsecond uncertainty and < 



1 km uncertainty in distance. For the instrumental delays, we used the known fiber delay values. 
The resulting delays have an accuracy of ~ 5 ns, a value smaller than the 10 ns sample time. The 
delay solution uncertainty is dominated instead by the delay range expected for sources distributed 



over the part of the lunar limb illuminated by the primary beam of the antennas (Fig. 10). The 
range of delay window widths varied with projected antenna separation, but was generally in the 
range from 80 ns to 140 ns. 



Gaussian statistics will generate many 'accidental' coincidences, as described in section 2.4, but 
the rate is greatly reduced as the threshold level and the number of antennas increases. For RESUN- 
A we chose a threshold value at = 4.23. With a ±140 ns delay window, a four-antenna accidental 
coincidence would occur once every 450 observing hours, i.e., an accidental trigger probability 
p ~ 0.1 in 45 hours. As a check of the pulse detection statistics, we sampled the antenna signals at 
a threshold value of at — 3.8 and analyzed a small time range using a large delay window (±200 



ns). With this delay width we expect 4-channel accidental coincidences every 10 minutes. Fig. 11 
shows an example of an accidental 4-antenna coincidence which disappeared when the correct delay 
windows were applied. 
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Illustration of Beam Location Longest Array Baseline (~ 5.3 km) 




0369 12 15 o 2 4 6 8 10 



Approximate RA from Lunar Center (arcmin) Time Since Start of Observation (hr) 



Fig. 10. — Maximum differential delays from lunar limb sources contained within the VLA primary 
beam for the longest baselines in Arrays 1 and 2 at all observing epochs. 



3.2.1. Pulse detection verification tests 

Since there are no natural sources of nanosecond pulses in the sky to verify the pulse detection 
scheme, we built a broad-band portable pulse generator and flew it 40 m above the center of the VLA 
using a helium-filled balloon. The portable unit generated 25 ns wide pulses which were amplified 
and radiated from an omni- directional antenna. We pointed three of the subarray antennas at the 
pulse generator and examined the delays between pulses using both realtime oscilloscope traces and 
the RESUN data acquisition system. (The fourth antenna could not be tested since the balloon 
payload was below the minimum antenna elevation of 8°). After applying the instrumental delays 
supplied by NRAO staff and calculating geometrical delays to the pulse generator, we verified that 
pulses were received at the correct delays within the expected delays to ±10 ns precision. 
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Fig. 11. — Example of four channel accidental trigger with a wide acceptance window. The > 3.8<r 
trigger signal was initiated in channel 0, while the triggers in channels 1-3 were within a ±200 ns 
tolerance window centered on the corresponding differential (geometric plus instrumental) delays 
to the Moon. 



4. RESUN-A Observations and Results 

We observed the east and west limbs of the Moon using two 4-antenna subarrays for a total 
of 50 hours in five 10 hour sessions between 11-20 Feb 2008. Each sub-array used four EVLA 
antennas with a 50 MHz bandwidth centered at 1450 MHz. Unfortunately, because of unexpected 
suppression of signals greater than 3a for two EVLA antennas in the east sub-array (found during 
post- processing) , this sub-array was deemed unusable and was dropped from the analysis. In 
addition, about five observing hours were lost as a result of bad weather, operator errors, and 
equipment malfunctions. 
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Fig. 12. — RESUN-A sky coverage for all five observing epochs. For each epoch, the dark green 
stripe is the Moon's path, while the shaded cyan region is the approximate effective aperture sky 
coverage. The instantaneous aperture solid angle is shown in the inset at bottom right. Background 
AGN are denoted by circles,with width inversely proportional to distance. 



4.1. Isotropic Flux Upper Limit 



We detected no 4-station coincident events greater than a 3.98cr threshold level during an total 
observation period of 44.95 hours. The corresponding upper limit to isotropic differential neutrino 
flux is 

(4.1) 



E dN/dE < 0.003 EeV knTV 1 sr 



at 90% confidence level in the neutrino energy range 21.6 < log[E(eV)] < 22.6. This is (marginally) 
lower than the lowest published upper limits for lunar detection experiments in this energy range, 



and in particular confirms the GLUE upper limit reported by Gorham et al. (2004). Fig. 13 shows 



the flux upper limit for RESUN-A vs neutrino energy, along with a number of previously reported 
upper limits from other lunar target searches. In addition, we show the projected flux sensitivity 
limits for RESUN-B and a future SKA-class array ( James Sz Protheroe|[20 09 ) with 1 km 2 effective 
area at an center observing frequency 200 MHz, bandwidth 100 MHz, and 200 hours observing time. 
Note that all upper limits have been calculated using the analytic effective aperture calculation 



described in Gayley et al. (2009). These may differ from upper limits reported using Monte Carlo 



calculations, i.e. the GLUE search ( Gorham et al.||2004 ) and nuMoon ( Scholten et al.||2009 ), where 



the authors report upper limits nearly lOx lower. This discrepancy is discussed in detail in Gayley 



et al. (2009) and James & Protheroe (2009) 
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Fig. 13 also shows expected neutrino fluxes for a variety of proposed isotropic neutrino source 
models. It is clear that the existing lunar searches have not yet reached sensitivity levels which 



could detect neutrinos from these models, but that RESUN-B will test Z-burst models (Fodor et al. 



2002) and that SKA will test topological defect models and may probe the high end of the GZK 



model predictions, depending on total observing time and effective aperture. 



4-1.1. Predicted Counts vs Model plots 



Fig. 14 illustrates an alternative scheme to plotting isotropic flux limits versus model predic- 
tions that takes into account the predicted model-dependent detection counts in a more quantitative 
manner. The figure plots the number of expected pulse detections in a fixed aperture and total 
observing time versus the pulse signal strength. The plot shows both the accidental rate for a given 
antenna array scheme and threshold (nearly vertical line) , and the expected counts for a variety of 
source models. The intersection between the accidental rate and the horizontal one count line is a 
fiducial point: A search that could detect neutrinos from a given model must have its fiducial point 
below the model. The cross-hatched area for each of the four example searches shown in the Figure 
is the region for which model counts exceed accidental counts. While RESUN-A does not have any 
cross-hatch areas (i.e, does not test any proposed isotropic models), RESUN-B will probe both halo 
(dash-dot) and extragalactic (bold dash-dot) Z-burst models, with ~1-10 expected counts during 
the 200 hour experiment. 



4.2. Point Source Flux Upper Limit 

Although RESUN-A did not have adequate sensitivity to test isotropic neutrino source models, 
the non-detection result can be used to determine upper limits for UHE neutrino point source (e.g. 



active galaxies) fluxes along the lunar path. Fig. 12 shows the sky coverage of RESUN-A over the 



entire five-day experiment, along with the sky positions of active galaxies (Veron-Cetty & Veron 



2006). For point sources along the lunar path shown in the figure, the upper limit to neutrino flux 



is (eqn. ^29) 



E 2 dN/dE < 3 x 10~ 4 EeV km _2 s _1 (4.2) 
at 90% confidence level in the neutrino energy range 21.6 < log[E(eV)] < 22.6. 



5. Summary 

In this paper, we present a complete, self-consistent description of the physical and technical 
considerations associated with searches for Cerenkov radio pulses from UHE neutrino interactions 
with the Moon. We describe the physics of the interaction, properties of the resulting Cerenkov 
radio pulse, detection statistics, effective aperture scaling laws, and derivation of upper limits for 
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Fig. 13. — UHE Neutrino isotropic flux upper limits determined by lunar neutrino searches: RESUN 
Phase -A (red line), and projected for RESUN Phase-B (blue line) along with previously reported 
upper limits (red dash-dot lines) from Parkes (Hankins et al. |1996| James et al. 2007), GLUE 
( Gorham et~al]|2004[ ) , Kalyazin ( |Beresnyak et al.||2005[ ), nuMoon ( |Scholten et~al"1|2009[ ). Note that 
all upper limits have been computed using the analytic aperture calculation of Gayley et al. (2009) 
and do not necessarily agree with Monte Carlo simulations. Also shown are cosmogenic neutrino 
flux predictions (James & Protheroe 2009 and refs.) from Z-bursts in the galactic halo and from 
extragalactic background (Fodor et al. 2002), GZK interactions with protons and with heavy ions 
(blue shaded band), topological defects (TD, green short dashes), and Waxman-Bahcall upper 



limits for no redshift evolution and with redshift evolution (WB, black short dashes, Waxman &; 



Bahcall (1999)). 
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Fig. 14. — Expected number of pulses exceeding trigger threshold versus detector sensitivity for 
four lunar neutrino detection experiments: (a) RESUN-A (45hr), (b) RESUN-B (200 hr), (c) 
nuMoon (20 hr), and (d) SKA (200 hr). In each case, the solid line is the accidental rate given 
by equation |2.15| while the labeled dotted and dashed lines are the expected rates for neutrino- 
production models [legend in panel (a)]. The shaded areas are regions where model counts exceed 
the accidental counts. In order to detect neutrinos from a given source model, the intersection 
between the accidental rate (nearly vertical line) and the unity count horizontal line must lie below 
the model's curve. 



isotropic and point source models. We report on initial results from the RESUN search, which 
uses multiple subarrays of the Expanded Very Large Array at 1.45 GHz pointing along the lunar 
limb. We detected no lunar pulses during 44.95 observing hours. The inferred isotropic neutrino 
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flux limit is comparable to the lowest published upper limits for lunar searches. The full RESUN 
search, with 200 hours observing time and an improved data acquisition scheme, will be be an order 
of magnitude more sensitive than previous lunar-target searches, and will test Z burst models of 
neutrino generation. 

We are grateful to Barry Clark, John Ralston, Peter Gorham, and Hallsie Reno for several 
illuminating discussions, Dan Mertley for technical assistance at the VLA, and the staff of CASPER 
for technical assistance with the data acquisition systems. The National Radio Astronomy Obser- 
vatory is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

REFERENCES 

Alvarez-Muhiz, J., Marques, E., Vazquez, R. A., & Zas, E. 2006, Phys. Rev. D, 74, 023007 

Alvarez-Muniz, J., James, C. W., Protheroe, R. J., & et al. 2008, International Cosmic Ray Con- 
ference, 5, 1507 

Alvarez-Muniz, J., Vazquez, R. A., & Zas, E. 2000, Physical Review D, 61 

Askaryan, G. A. 1962, Journal of the Physical Society of Japan Supplement, 17, C257+ 

Barwick, S. W., et al. 2006, Physical Review Letters, 96, 171101 

Beresnyak, A. R., Dagkesamanskii, R. D., Zheleznykh, I. M., Kovalenko, A. V., & Oreshko, V. V. 
2005, Astronomy Reports, 49, 127 

Butler, B. J., McKinnon, M. M., Perley, R. A., & Dewdney, P. E. 2006, IAU Special Session, 1 

Cronin, J. W., Gaisser, T. K., & Swordy, S. P. 1997, Scientific American, 276, 32 

Fodor, Z., Katz, S. D., & Ringwald, A. 2002, Physical Review Letters, 88 

Fodor, Z., Katz, S. D., & Ringwald, A. 2002, ArXiv High Energy Physics - Phenomenology e-prints 
Gayley, K. C, Mutel, R. L., & Jaeger, T. R. 2009, ArXiv e-prints, 0904 
Gorham, P. W., et al. 2009, Physical Review Letters, 103 
— . 2007, Physical Review Letters, 99 

Gorham, P. W., Hebert, C. L., Liewer, K. M., Naudet, C. J., Saltzberg, D., & Williams, D. 2004, 
Physical Review Letters, 93, 041101 

Gorham, P. W., Saltzberg, D. P., Schoessow, P., Gai, W., Power, J. G., Konecny, R., & Conde, 
M. E. 2000, Phys. Rev. E, 62, 8590 



- 25 - 



Greisen, K. 1966, Physical Review Letters, 16, 748 
Gusev, G. A., et al. 2006, Cosmic Research, 44, 19 

Hankins, T. H., Ekers, R. D., & O'Sullivan, J. D. 1996, MNRAS, 283, 1027 

James, C. W., Crocker, R. M., Ekers, R. D., Hankins, T. H., O'Sullivan, J. D., k Protheroe, R. J. 
2007, mnras, 379, 1037 

James, C. W., Ekers, R. D., Alvarez-Muiz, J., Protheroe, R. J., McFadden, R. A., Phillips, C. J., 
& Roberts, P. 2009, Nuclear Instruments and Methods in Physics Research A, 604 

James, C. W., & Protheroe, R. J. 2009, Astroparticle Physics, 31, 392 

James, C. W., & Protheroe, R. J. 2009, Astroparticle Physics, 30, 318 

Kravchenko, I. e. a. 2006, Phys. Rev. D, 73, 082002 

Lehtinen, N. G., Gorham, P. W., Jacobson, A. R., & Roussel-Dupre, R. A. 2004, Phys. Rev. D, 69, 
013008 

Meeus, J. 1998, Astronomical algorithms (Richmond, VA: Willmann-Bell) 

Miocinovic, P., Field, R. C, Gorham, P. W., Guillian, E., Milincic, R., Saltzberg, D., Walz, D., & 
Williams, D. 2006, Physical Review D, 74 

Olhoeft, G. R., & Strangway, D. W. 1975, Earth and Planetary Science Letters, 24, 394 

Panda, S., Mohanty, S., Padmanabhan, J., & Stal, O. 2007, Journal of Cosmology and Astro- 
Particle Physics, 11, 22 

Parsons, A., et al. 2008, Publications of the Astronomical Society of the Pacific, 120, 1207 

Reno, M. H. 2005, Nuclear Physics B Proceedings Supplements, 143, 407 

Saltzberg, D., et al. 2001, Physical Review Letters, 86, 2802 

Scholten, O., Bacelar, J., Braun, R., de Bruyn, A. G., Falcke, H., Stappers, B., & Strom, R. G. 
2006, Astroparticle Physics, 26, 219 

Scholten, O., et al. 2009, Nuclear Instruments and Methods in Physics Research A, 604, 102 

Scholten, O., et al. 2009, ArXiv e-prints, 0910.4745S 

Shepard, M. K., Brackett, R. A., & Arvidson, R. E. 1995, J. Geophys. Res., 100, 11709 
Stanev, T. 2004, Physics Letters B, 595, 50 

Veron-Cetty, M. P., & Veron, P. 2006, Astronomy and Astrophysics, 455, 773 



- 26 - 



Waxman, E. 2009, Nuclear Physics A, 827, 15 

Waxman, E., & Bahcall, J. 1999, Physical Review D, 59 

Williams, D. R. 2004, PhD thesis, University of California, Los Angeles 

Yamamoto, T. 2008, International Cosmic Ray Conference, 4, 335 



This preprint was prepared with the AAS IATjgX macros v5.2. 



